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Introduction and Theory 


The desirable properties of ceramics at high temperatures have generated interest in their use for 
structural applications such as in advanced turbine systems. Design lives for such systems can 
exceed 10,000 hours. Such long life requirements necessitate subjecting the components to 
relatively low stresses. The combination of high temperatures and low stresses typically places 
failure for monolithic ceramics in the creep regime (ref 1). 

The objective of this work is to present a design methodology for predicting the lifetimes of 
structural components subjected to multiaxial creep loading. This methodology utilizes 
commercially available finite element packages and takes into account the time varying creep 
stress distributions (stress relaxation). In this methodology, the creep life of a component is 
divided into short time steps, during which, the stress and strain distributions are assumed 
constant. The damage, D, is calculated for each time step based on a modified Monkman-Grant 
creep rupture criterion (ref 2). For components subjected to predominantly tensile loading, failure 
is assumed to occur when the normalized accumulated damage at any point in the component is 
greater than or equal to unity. 

Some ceramic components, such as vanes and rotors, are subjected to concurrent tensile and 
compressive stress fields. For such components, failure generally starts at or near the most highly 
stressed point and subsequently propagates across the section. The creep rupture life for members 
subjected to concurrent tensile and compressive loading is divided into two stages. The first is 
called the stage of latent failure (damage initiation). During this stage, the damage accumulates 
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until it reaches unity at some point within the component, and failure begins. Hence, this portion 
of life for members subjected to concurrent tensile and compressive loading represents the entire- 
predicted life for predominantly tensile components. Damage due to compressive stresses is 
assumed to be negligible, although in this methodology, it can be accounted for very easily if it is 
determined to be of any significance. Subsequently, the second stage, named the damage 
propagation takes place. During this stage a damage front defined by the condition D=1 will 
travel through the body or surface of the component. Component failure occurs at the end of this 
stage when its total load carrying capacity is expended. In CARESICreep (Ceramics Analysis 
and Reliability Evaluation of Structures/Creep), this means that failure is assumed to occur when 
D=1 at the periphery of the expanded critical damage zone. The corresponding time will be the 
creep rupture life for that component. This size of the critical damage zone corresponding to 
creep rupture failure varies depending on the load conditions and component configuration. One 
estimation for the duration of the propagation stage is to assume it equal to the time it takes the 
damage zone to penetrate the initial tensile stressed portion of the structure. 

Two examples were chosen to demonstrate the viability of the creep life prediction methodology 
presented above. The integrated design code CARES, I Creep (Fig. 7) which utilizes this damage 
accumulation model was used for this purpose. The first example, silicon nitride NCX-5100 
notched tensile specimen, represents the application of this approach to predicting the creep 
rupture life for components subjected to multiaxial predominantly tensile creep loading. The 
second example, siliconized silicon carbide KX01 bend specimen, represents the application of 
this approach to components subjected to uniaxial tensile-compressive creep loading. It was 
found that the methodology described in this paper yielded good creep rupture life predictions for 
both examples, given the amount of scatter usually found in the creep rupture life data. 

Examples 

D Notched Tensile Specimen: 

The creep experiments were conducted on two types of silicon nitride NCX-5100 specimens. 
First, smooth tensile tests were investigated in order to characterize the creep response of the 
silicon nitride. Second, experiments on notched tensile bars (Fig. 8) provided a multiaxial stress 
state, especially near the notch root, where the creep life may be predicted from data obtained 
from the smooth tensile specimens. Creep data used in this example were obtained from ref 3. 

Using CARESICreep, the secondary creep rate parameters, C 7 , C s , and C I0 , were determined to 
be 7.858xl0' 26 /Pa 6 75 hour, 6.75, and 127560°K, respectively. Figure 9 shows the experimental 
secondary creep strain rate vs. the analytical secondary creep strain rate. This figure shows that 
the data scatters relatively close to the 45° line, indicating that the Norton secondary creep rate 
model (Fig. 9) was successful in describing the secondary creep behavior of the NCX-5 100 
material. 

Four notched bars were tested (ref 3). The reduced section average stresses were 105, 120, 135, 
and 150 MPa. The maximum principal stress distribution for the 120 MPa reduced average 
section stress specimen as a fimction of time is shown in Fig. 10. A multiaxial stress state exists 
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in the vicinity of the notch root, while the stress state away from the notch is constant and 
uniaxial. When the load is initially applied, time=0, the maximum principal stress is located at 
the root of the notch. As time progresses, the stress relaxes and, the local maximum moves into 
the interior of the notched bar. This stress relaxation, which occurs due to nonlinear creep 
deformation, will influence the damage calculations. This is because it was found that the 
location of the maximum cumulative damage also moves away from the surface as time elapses, 
indicating that failure could originate at the interior of the specimen. 

The Monkman-Grant criterion in association with the damage model was used to predict the 
creep life for the notched specimens. Since this component is entirely subjected to tensile stress 
state, failure was assumed to occur when D first reached unity at any point. The Table shown in 
Fig. 1 1 gives a summary of the failure predictions for the notched bars as a function of their 
reduced section stress. The predicted lives using the methodology described in this paper and 
computed via CAKES/Creep for these specimens compare well to the experimental failure times. 
Figure 1 1 also shows a cumulative damage map for the 120 MPa specimen after 80 hours. The 
cumulative damage is equal to one and is located near to but not directly the root of the notch. 


21 Bend bar 

This example is presented to demonstrate the methodology discussed in this paper for predicting 
the creep life of ceramic components subjected to simultaneous tensile/compressive stresses. 
Wiederhom et al. (ref 4) conducted creep testing on the KX01 siliconized silicon carbide material 
at 1300°C using flexure, tensile and compressive specimens. They found that this material spent 
most of its life in the secondary creep region. Further, they found that the KX01 material displays 
significant asymmetric creep behavior, and that in both tension and compression, the creep rate 
displayed a bilinear power (Norton) law behavior with a transition point at a threshold stress 
(Fig. 13). In this example, CARES/Creep was used to predict the creep lives of four point bend 
specimens using the tensile and compressive creep and creep rupture properties of the material. 
These predicted lives were then compared to the experimental ones. 

Figure 14 displays the evolution of stress as function of time (stress relaxation) for the bar tested 
at 250 MPa. The stress distribution in the bend bar converged to its stationary value ten hours 
after the load was applied. This figure also shows how the neutral axis shifts towards the 
compressive region as the specimen creeps. 

Figure 15 show the evolution of damage in the flexure bar tested at 250 MPa initially, at t=14.5 
hours corresponding to the time when the damage first reached unity (latent stage of failure), and 
at t=80 hours corresponding to the time when the final predicted failure occurred (end of damage 
propagation stage). In this analysis, it was assumed that failure would occur when the initial 
portion of the bend bar stressed in tension is damage (half the depth). The table shown in Fig. 16 
indicates that excellent agreement between experimental and predicted creep lives exists for the 
bend bars. 

Several advantages are apparent to this creep rupture life methodology. First, this methodology 
yields a cumulative damage map for the component showing the critical locations where failure 
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would originate. This capability is very helpful in the redesign of such components. In creep type 
loading applications, it is not a trivial task to predict the location of failure since the multiaxial 
stress components redistribute as time elapses. Thus, failure will not necessarily occur at the 
location where stresses are highest at the beginning of loading or at the time of failure, but can 
take place elsewhere, as shown by the notched tensile specimen. Second, this design 
methodology is capable of incorporating the primary creep strain effect into the analysis 
(influences the stress state) which could predict shorter lives (conservative predictions) compared 
to when only the secondary creep stain effect is used. Third, any equivalent stress criterion can 
be used to predict the component’s life. As multiaxial creep data emerge in the future and we 
understand better how ceramics fail under such applications, the CAKES! Creep code can be 
modified accordingly. Fourth, any creep rupture criterion (Continuum damage mechanics, 
Larson-Miller, minimum commitment method, etc.) can be utilized to compute the damage and 
predict life. 
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Objective 


To develop an integrated design program for predicting the lifetime of 
structural ceramic components subjected to multiaxial creep loads. 
This methodology utilizes commercially available finite element 
packages and takes into account the time varying creep stress 
distribution (stress relaxation). 


Fig. 1 


Constitutive Creep Laws 


Total Creep Strain 

e cr = e Primary + Secondary 

Primary Creep Strain Increment — Bailey-Norton Law 
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Secondary Creep Strain Increment — Norton Law 



where cr is the stress 
t is the time 

T is the absolute temperature 

Cj are experimentally determined constants 
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Multiaxial Creep Model 


Flow Rule 
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where Sy is the deviatoric stress 


Fig. 3 


Effective Stress 

CT e = ^={(oi-o 2 ) 2 + ( a 2 - o 3 ) 2 + ( o 3 - cj -|) 2 } 1/2 
Effective Secondary Creep Strain 



where o-j , o 2 . and o 3 are the principal stresses 


Fig. 4 
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Creep Rupture 


Monkman-Grant 


tf = b-| £ s 


-b 2 


Modified Monkman-Grant 

Temperature Stratification 


In t^ = d-j — d g In Eg + 



where bj and dj are constants 


Fig. 5 


Damage 


Damage Assessment 

0 < D < 1 


where D = 0 for an undamaged component 
D = 1 for a failed component 


N 

D = E 


At: 


j=1 exp [ d 1 - d 2 In e s . + d 3 /Tj ] 


where At: is the duration of the i th time step 
In is the number of time steps 


Fig . 6 
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Integrated Design Program 



Fig. 7 

CARES /Creep Multiaxial Tensile Creep Benchmark 

Notched Tensile Specimen - Norton Company 

Axisymmetric Model 
1047 PLANE82 elements 


NCX-5100 Silicon Nitride 
Axisymmetric Model 
Temperature = 1370°C 
Constant Applied Load 


Height = 8.9 cm 
Model Height = 0.72 cm 
Diameter = 0.32 cm 
Notch Radius = 0.032 cm 



Fig. 6 
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Parameter Estimation 

Silicon Nitride - N CX-5100 


37 Tensile Specimen s 
Temperatures = 1275-1425 °C 
Applied Load = 120-250 MPa 
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Analytical creep rate (1/hr) 


Fig. 9 


Stress Relaxation Near the Notch Root 

Nominal Reduced Section Stress=120 MPa 
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Notched Tensile Specimen 


Damage Map 



Results 
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Fig. 11 


CARES /Creep Uniaxial Tensile - Compressive 
Creep Benchmark 

Flexural Specimen - NIST 


KX01 Siliconized Silicon Carbide 
Plane Stress Model 
Temperature = 1300°C 
Constant Applied Load 



Height = 4.0 mm 
Length = 50 mm 
Depth = 3 mm 


Fig. 12 
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Creep rate (1/hr) 


Parameter Estimation 

Siliconized Silicon Carbide - KX01 
Uniaxial Tension and Compression Tests at 1300°C 



Material displays asymmetric bilinear creep behavior in tension and compression 

Fig. 13 

Stress Relaxation 

Maximum Elastic Stress=250 MPa 



Fig. 14 
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Maximum Elastic Stress=250 MPa 



MPa 


Time = 14.5 hrs 



.113095 

.226189 

.339284 

.452379 

.565473 

.678568 

.791663 

.904757 

1.018 



Time = 80 hrs 


□ 


0 

1 

>1 


Fig. 15 


Flexure Specimen 


Results 
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Summary and Conclusions 


A general purpose creep life prediction code, CARES/ Creep, has been 
developed. It is integrated with ANSYS finite element software, and can be 
used to design monolithic ceramic components. 

Creep life was predicted based on accumulated damage and the 
Monkman-Grant and the Modified Monkman-Grant failure criteria. 


The creep life was well predicted for components subjected to multiaxial 
tensile and uniaxial tensile-compressive stress states. 

Life prediction for a component with simultaneous tensile-compressive 
stress states is a two step process involving damage initiation and 
propagation. 

This creep life prediction methodology results in damage maps showing 
critical regions and, hence, can be used in the components design. 


Fig. 17 


Focus of Future Research 


Creep Response 

Multiaxial Tensile-Compressive Modeling 
Continuum Damage Mechanics 

Creep Rupture 

Probabilistic Methods 
Creep-Fatigue Interaction 
Integration with CARES/ Life 


Fig. 18 
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